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In this paper we investigated structural, electronic and magnetic properties of 3d (light) transition 
metal (TM) atomic chains using first-principles pseudopotential plane wave calculations. Periodic 
linear, dimerized linear and planar zigzag chain structures and their short segments consisting of 
finite number of atoms have been considered. Like Cu, the periodic, linear chains of Mn, Co and Ni 
correspond to a local shallow minimum. However, for most of the infinite periodic chains, neither 
linear nor dimerized linear structures are favored; to lower their energy the chains undergo a struc- 
tural transformation to form planar zigzag and dimerized zigzag geometry. Dimerization in both 
infinite and finite chains are much stronger than the usual Peierls distortion and appear to depend on 
the number of 3d-electrons. As a result of dimerization, a significant energy lowering occurs which, 
in turn, influences the stability and physical properties. Metallic linear chain of Vanadium becomes 
half-metallic upon dimerization. Infinite linear chain of Scandium also becomes half-metallic upon 
transformation to zigzag structure. An interplay between the magnetic ground state and atomic 
as well as electronic structure of the chain has been revealed. The end effects influence the geom- 
etry, energetics and magnetic ground state of the finite chains. Structure optimization performed 
using noncollinear approximation indicates significant differences from the collinear approximation. 
Variation of the cohesive energy of infinite and finite-size chains with respect to the number of 
Sri-electrons are found to mimic the bulk behavior pointed out by Friedel. The spin-orbit coupling 
of finite chains are found to be negligibly small. 

PACS numbers: 73.63.Nm, 75.50.Xx, 75.75.+a 



I. INTRODUCTION 

Fabrication of nanoscale structures, such as quan- 
tum dots, nanowires, atomic chains and functionalized 
molecules have made a great impact in various fields 
of science and technology^! 3 -^ Size and dimensionality 
have been shown to strongly affect the physical and chem- 
ical properties of matter^ Electrons in lower dimension- 
ality undergo a quantization which is different from that 
in the bulk materials.— i 7 -^ In nanoscale size, the quantum 
effects, in particular the discrete nature of electronic en- 
ergies with significant level spacing are pronounced. 

Suspended monoatomic chains being an ultimate one- 
dimensional (ID) nanowire have been produced and their 
fundcmcntal properties both theoretically and experi- 
mentally investigated^^ 2 .^^^^^ Ballistic 
electron transport^ with quantized conductance at room 
temperature has been observed in metallic nanowires 
Moreover, magnetic and transport properties become 
strongly dependent on the details of atomic configura- 
tion. Depending on the type and position of a foreign 
atom or molecule that is adsorbed to a nansostructure, 
dramatic changes can occur in the physical properties. 3 
Some experimental studies, however, aimed at producing 
atomic chains on a substrate^ Here the substrate-chain 
interaction can enter as a new degrees of freedom to in- 
fluence the physical properties. 

Unlike the metal and semiconductor chains, not 
many theoretical studies are performed on transition 
meta l 21 i 22 ' 23 i 24 (TM) monatomic chains. TM monatomic 



chains have the ability to magnetize much more than the 
bulk. 25 Large exchange interactions of TM atoms in the 
bulk are overcame by the large electron kinetic energies, 
which result in a nonmagnetic ground state with large 
bandwidth. On the other hand, geometries which are 
nonmagnetic in bulk may have magnetic ground states 
in monatomic chains^ In addition, it is predicted that 
the quantum confinement of electrons in metallic chains 
should result in a magnetic ground state and also a super 
paramagnetic state for some of the TM chains— at finite 
temperatures. The central issue here is the stability of 
the chain and the interplay between ID geometry and 
magnetic ground state i 21 i 24 

From the technological point of view, TM monatomic 
chains are important in spin dependent electronics, 
namely spintronics. 27 While most of the conventional 
electronics is based on the transport of information 
through charges, future generation spintronic devices will 
take the advantage of the electron spin to double the 
capacity of electronics. It has been revealed that TM 
atomic chains either suspended or adsorbed on a ID sub- 
strate, such as carbon nanotubes or Si nanowires can 
exhibit high spin-polarity or half-metallic behavior rele- 
vant the for spin- valve effect 4 Recently, first-principles 
pseudopotential calculations have predicted that the fi- 
nite size segments of linear carbon chains capped by spe- 
cific 3d-TM atoms display an interesting even-odd dis- 
parity depending on the number of carbon atoms. For 
example, CoC n Co linear chain has an antiferromagnetic 
ground state for even n, but the ground state changes to 
ferromagnetic for odd n. Even more interesting is the fer- 
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romagnetic excited state of an antiferromagnetic ground 
state can operate as a spin-valve when CoC„Co chain is 
connected to metallic electrodes from both ends.— 

As the chain length decreases, finite-size effects dom- 
inate the magnetic and electronic properties . 21 ' 29 When 
compared with the infinite case, these are less stable to 
thermal fluctuations^ Additional effects on the behavior 
of nanoparticle are their intrinsic properties and the in- 
teraction between them i 29 i 30 i 31 i 32 Effects of noncollinear 
magnetism have to be taken into account as wel L 33 ' 34 ' 35 
The end atoms also exhibit different behavior with re- 
spect to the atoms close to the middle of the structured 

In this paper, we consider infinite, periodic chains of 
3d-TM atoms having linear and planar zigzag structures 
and their short segments consisting of finite number of 
atoms. For the sake of comparison, Cu and Zn chains 
also included in our study. All the chain structures dis- 
cussed in this paper do not correspond to a global mini- 
mum, but may belong to a local minimum. The infinite 
and periodic geometry is of academic interest and can 
also be representative for very long monatomic chains. 
The main interest is, however, the short segments com- 
prising finite number of TM atoms. We examined the 
variation of energy as a function of the lattice constant 
in different magnetic states, and determined the stable 
infinite and also finite-size chain structures. We inves- 
tigated the electronic and magnetic properties of these 
structures. Present study revealed a number of proper- 
ties of fundamental and technological interest: The linear 
geometry of the infinite, periodic chain is not stable for 
most of the 3d-TM atoms. Even in the linear geometry, 
atoms are dimerized to lower the energy of the chain. We 
found that infinite linear Vanadium chains are metallic, 
but become half-metallic upon dimerization. The pla- 
nar zigzag chains are more energetic and correspond to 
a local minimum. For specific TM chains, the energy 
can further be lowered through dimer formation within 
planar zigzag geometry. Dramatic changes in electronic 
properties occur as a result of dimerization. Magnetic 
properties of short monatomic chains have been investi- 
gated using both collinear and noncollinear approxima- 
tion, which are resulted in different net magnetic moment 
for specific chains. Spin-orbit coupling which are calcu- 
lated for different direction of applied magnetic field have 
been found to be negligibly small. 



structures has been chosen as 0.1 eV for geometry 
relaxations and as 0.01 eV for accurate energy band 
and density of state (DOS) calculations. As for finite 
structures, the width of smearing is taken as 0.01 eV. 
We treated the chain structures by supercell geometry 
(with lattice parameters a sc , b sc , and c sc ) using the 
periodic boundary conditions. A large spacing (~ lOA) 
between adjacent chains has been assured to prevent 
interactions between them. In single cell calculations of 
infinite systems, c sc has been taken to be equal to the 
lattice constant of the chain. The number of plane waves 
used in expanding Bloch functions and k-points used in 
sampling the Brillouin zone (BZ) have been determined 
by a series of convergence tests. Accordingly, in the 
self-consistent potential and total energy calculations the 
BZ has been sampled by (1x1x41) mesh points in k-space 
within Monkhorst-Pack scheme.— A plane-wave basis 
set with kinetic energy cutoff ?i 2 |k + G\ 2 /2m = 350 eV 
has been used. In calculations involving PAW potentials, 
kinetic energy cutoff is taken as 400 eV. All atomic 
positions and lattice constants (c sc ) have been optimized 
by using the conjugate gradient method where total en- 
ergy and atomic forces are minimized. The convergence 
is achieved when the difference of the total energies of 
last two consecutive steps is less than 10 -5 eV and the 
maximum force allowed on each atom is 0.05 eV/A. 
As for finite structures, supercell has been constructed 
in order to assure ~ lOA distance between the atoms 
of adjacent finite chain in all directions and BZ is 
sampled only at the T-point. The other parameters of 
the calculations are kept the same. The total energy of 
optimized structure (£r) relative to free atom energies 
is negative, if they are in a binding state. As a rule, 
the structure becomes more energetic (or stable) as its 
total energy is lowered. Figure [1] describes various chain 
structures of TM atoms treated in this study. These are 
infinite periodic chains and segments of a small number 
of atoms forming a string or a planar zigzag geometry. 
The stability of structure-optimized finite chains are 
further tested by displacing atoms from their equilibrium 
positions in the plane and subsequently reoptimizing 
the structure. Finite-size clusters of TM atoms are kept 
beyond the scope of this paper. 



II. METHODOLOGY 

We have performed first-principles plane wave 
calculation !] 37 ' 38 within Density-Functional Theory 
(DFT) 39 using ultra-soft pseudopotentials4^ We also 
used PAW— potentials for noncollinear and non- 
collinear spin-orbit calculations of finite chains. The 
exchange-correlation potential has been approximated 
by generalized gradient approximation (GGA)42 For 
partial occupancies, we have used the Methfessel-Paxton 
smearing method.— The width of smearing for infinite 



III. INFINITE AND PERIODIC CHAIN 
STRUCTURES 

In Fig. [2 we present the energy versus lattice con- 
stant of various infinite and periodic chain structures 
(described in Fig. [T]) in different magnetic states. In cal- 
culating the ferromagnetic (FM) state, the structure is 
optimized each time using a spin-polarized GGA calcula- 
tions starting with a different preset magnetic moment in 
agreement with Hund's rule. The relaxed magnetic mo- 
ment yielding to the lowest total energy has been taken 
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FIG. 1: (Color online) Various structures of 3d- TM atomic 
chains, (a) Infinite and periodic structures; L: The infinite 
linear monatomic chain of TM atom with lattice constant, 
c. LD: The dimerized linear monatomic chain with two TM 
atoms in the cell, e is the displacement of the second atom 
from the middle of the unit cell. ZZ: The planar zigzag 
monatomic chain with lattice parameter c and unit cell hav- 
ing two TM atoms, ci ~ c 2 and 59° < a < 62°. ZZD: 
The dimerized zigzag structure ci 7^ C2. WZ: The wide angle 
zigzag structure ci ~ C2, but a > 100°. (b) Various chain 
structures of small segments consisting of finite number (n) 
of TM atoms, denoted by (TM)„. 



as the FM state of the chain. For the antiferromagnetic 
(AFM) state, we assigned different initial spins of oppo- 
site directions to adjacent atoms and relaxed the struc- 
ture. We performed spin-unpolarized GGA calculations 
for the nonmagnetic (NM) state. The energy per unit 
cell relative to the constituent free atoms is calculated 
from the expression, E = [NE a — Et], in terms of the 
total energy per unit cell of the given chain structure 
for a given magnetic state (Et) and the ground state 
energy of the free constituent TM atom, E a . N is the 
number of TM atom in the unit cell, that is N=l for L, 
but N=2 for LD, ZZ, and ZZD structures. The mini- 
mum of E is the binding energy. By convention E^ < 
corresponds to a binding structure, but not necessary 
to a stable structure. The cohesive energy per atom is 
E c = —Eb/N. Light transition metal atoms can have dif- 



ferent structural and magnetic states depending on the 
number of their 3d electrons. For example, Sc having a 
single 3d electron, has a shallow minimum corresponding 
to a dimerized linear chain structure in the FM state. If 
the L structure is dimerized to make a LD structure, the 
energy of the chain is slightly lowered. Other linear struc- 
tures, such as linear NM, and AFM have higher energy. 
More stable structure ZZ is, however, in the FM state. 
The situation is rather different for Cr, Fe, and Mn. For 
example, Cr have LD and more energetic ZZD structures 
in the AFM state. It should be noted that in the dimer- 
ized linear chain structure of Cr the displacement of the 
second atom from the middle of the unit cell, e, is rather 
large. Apparently, the dimerization is stronger than the 
usual Peierls distortion. As a result, the nearest neigh- 
bor distance, (c — e), is much smaller than the second 
nearest neighbor distance, (c + e). This situation poses 
the question whether the interaction between the adja- 
cent dimers are strong enough to maintain the coherence 
of the chain structure. We address this question by com- 
paring the energies of individual dimers with the chain 
structure. The formation of the LD structure is ener- 
getically more favorable with respect to individual dimer 
by 0.54 eV per atom. Furthermore, the charge density 
contour plots presented in Fig. [3] indicate a significant 
bonding between adjacent dimers. Nevertheless, the LD 
structure has to transform to more energetic ZZD struc- 
ture. The zigzag structures in the AFM, FM and NM 
states have minima at higher binding energies and hence 
are unstable. The linear and LD Fe chains have a local 
minimum in FM state. More stable ZZ and ZZD struc- 
tures in FM state have almost identical minima in lower 
binding energy. The most stable chain structure of Mn 
among ones described in Fig. [2] is ZZ in FM state. It 
is also saliency to note that both Fe and Mn chains in 
NM state undergo a structural transformation from ZZ to 
WZ structure. It is noted from Fig. [2] that the structure 
of 3d-TM atomic chains are strongly dependent on their 
magnetic state. Optimized structural parameters, cohe- 
sive energy, magnetic state and net magnetic moment 
of infinite linear and zigzag structures are presented in 
Table U and Table [Til respectively. 



In Fig. |4] and Fig. [5] we compare the nearest neighbor 
distance and the average cohesive energy of the linear 
and zigzag chain structures with those of the bulk met- 
als and plot their variations with respect to their number 
of 3d electrons of the TM atom. The nearest neighbor 
distance in the linear and zigzag structures are smaller 
than that of the corresponding bulk structure, but dis- 
play the similar trend. Namely, it is large for Sc having 
a single 3d electron and decreases as the number of 3d 
electrons, i.e. N4, increases to four. Mn is an exception, 
since the bulk and the chain structure show opposite be- 
havior. While the nearest neighbor distance of bulk Mn 
is a minimum, it attains a maximum value in the chain 
structure. Owing to their smaller coordination number, 
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TABLE II: The calculated values for the planar zigzag struc- 
tures (ZZ and ZZD): The lattice constant, c (in A); the first 
nearest neighbor, ci (in A); the second nearest neighbor, C2 (in 
A); angle between them, a (in degrees); the cohesive energy, 
E c (in eV/atom); the magnetic ground state, MGS; and the 
total magnetic moment, fi per unit cell (in Bohr magnetons, 
fis) obtained within collinear approximation. 



FIG. 2: (Color online) The energy versus lattice constant, c, 
of various chain structures in different magnetic states. FM: 
ferromagnetic; AFM: antiferromagnetic; NM: nonmagnetic; 
FMD: ferromagnetic state in the linear or zigzag dimerized 
structure; AFMD: antiferromagnetic state in the dimerized 
linear or zigzag structure. The energy is taken as the energy 
per unit cell relative to the constituent free atom energies 
in their ground state (See text for definition). In order to 
compare the energy of the L structure with that of the LD, the 
unit cell (and also lattice constant) of the former is doubled 
in the plot. Types of structures identified as L, LD, ZZ, ZZD, 
WZ are describes in Fig. [T] 





Sc 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


c 


6.0 


4.9 


4.5 


4.4 


2.6 


4.6 


2.1 


2.2 


2.3 


2.6 


e 


0.38 


0.52 


0.51 


0.66 


0.0 


0.21 


0.0 


0.0 


0.0 


0.0 


E c 


1.20 


1.83 


1.86 


1.40 


0.76 


1.81 


2.10 


1.99 


1.54 


0.15 


MGS 


FM 


FM 


FM 


AFM 


AFM 


FM 


FM 


FM 


NM 


NM 




1.74 


0.45 


1.00 


±1.95 


±4.40 


3.32 


2.18 


1.14 


0.0 


0.0 



TABLE I: The calculated values for linear structures (L and 
LD). The lattice constant, c (in A); the displacement of the 
second atom in the unit cell of dimerized linear structure, e (in 
A); cohesive energy, E c (in eV/atom); the magnetic ground 
state, MGS; and the total magnetic moment, /j, per unit cell 
(in Bohr magnetons, ps) obtained within collinear approxi- 
mation. 
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FIG. 3: (Color online) Charge density contour plots for the 
linear (L) and the dimerized linear structure (LD) of Cr 
monatomic chains. Interatomic distances are indicated. Con- 
tour spacings are equal to Ap — 0.0827e/A 3 . The outermost 
contour corresponds to Ap = 0.0827e/A 3 . 



minima in the Born-Oppenheimer surface. 



chain structures have smaller cohesive energy as com- 
pared to the bulk crystals as shown in Fig. [5] However, 
both L (or LD if it has a lower energy) and ZZ (or ZZD 
if it has a lower energy) also show the well-known dou- 
ble hump behavior which is characteristics of the bulk 
TM crystals. Earlier, this behavior was explained for the 
bulk TM crystals4^^ The cohesive energy of zigzag 
structures are generally ~ 0.7 eV larger than that of the 
linear structures. However, they are 1-2 eV smaller than 
that of the bulk crystal. This implies that stable chain 
structures treated in this study correspond only to local 



We note that spin-polarized calculations are carried 
out under collinear approximation. It is observed that 
all chain structures presented in Table U and Table |TT] 
have magnetic state if Nd < 9. Only Cr and Mn 
linear chain structures and Cr zigzag chain structure 
have an AFM lowest energy state. The binding energy 
difference between the AFM state and the FM state, 
AE = E^ FM - E[ M , is calculated for all light TMs. 
Variation of AE with Nd is plotted in Fig. [51 We see 
that only Cr ZZ and ZZD chains have an AFM lowest 
energy state. AE of Fe is positive and has the largest 
value among all 3d-TM zigzag chains. Note that AE 
increases significantly as a result of dimerization. 
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FIG. 4: (Color online) Variation of the nearest neighbor dis- 
tance of 3d-TM atomic chains and the bulk structures. For 
the linear and zigzag structures the lowest energy configura- 
tion (i.e. symmetric or dimerized one) has been taken into 
account. Experimental values of the bulk nearest neighbor 
distances have been taken from Ref. [431. 




Sc Ti V Cr Mn Fe Co Ni Cu Zn 

FIG. 5: (Color online) Variation of the cohesive energy, E c 
(per atom), of 3d-TM monatomic chains in their lowest en- 
ergy linear, zigzag and bulk structures. For the linear and 
zigzag structures the highest cohesive energy configuration 
(i.e. symmetric or dimerized one) has been taken into ac- 
count. Experimental values of the bulk cohesive energies have 
been taken from Ref. 14711 




<|-D.2 



Sc Ti V Cr Mn Fe Co Ni Cu Zn 



FIG. 6: (Color online) Variation of the binding energy dif- 
ference, AE (per atom) between the lowest antiferromagnetic 
and ferromagnetic states of 3d-TM monatomic chains. Open 
squares and filled circles are for the symmetric zigzag ZZ and 
dimerized zigzag ZZD chains respectively. 



are Mn, Cr, and V chains. The linear Mn chain has an 
AFM state, where spin-up (majority) and spin-down (mi- 
nority) bands coincide. Chromium L and LD structures 
are AFM semiconductors. Vanadium is a ferromagnetic 
metal for both spins, but becomes half- metallic upon 
dimerization. In half-metallic state, the chain has integer 
number of net spin in the unit cell. Accordingly, Vana- 
dium chain in the LD structure is metallic for one spin 
direction, but semiconducting for the other spin direc- 
tion. Hence, the spin polarization at the Fermi level, i.e. 
P = [Df (Ep) - D l (E F )]/[D t (E F ) + D L {E F )} is 100%. 
Bands of Cu and Zn with filled 5c?-shell in nonmagnetic 
state are in agreement with previous calculations.— In 
Fig. [5J the chain structures in the first column have only 
ZZ structure. The zigzag chains in the second column 
are transformed to a lower energy (i.e. more energetic) 
ZZD structure in the third column. The ZZ chain of Sc 
is stable in a local minimum and displays a half-metallic 
character with 100% spin-polarization at the Fermi level. 
Accordingly, a long segment of ZZ chain of Sc can be 
used as a spin-valve. Ti, Mn, and Ni in their stable 
zigzag structures are FM metals. The stable ZZD struc- 
ture of Fe and Co chains are also FM metals. The ZZ 
and relatively lower energy ZZD structure of V chain are 
nonmagnetic. Both ZZ and ZZD structures of Cr are in 
the AFM state. 



Having discussed the atomic structure of 5d-TM 
chains, we next examine their electronic band structure. 
In Fig. the chain structures in the first column do 
not dimerize. The linear chains placed in the third col- 
umn are dimerized and changed from the L structure 
in the second column to form the LD structure. Most 
of the linear structures in Fig. [7] display a FM metallic 
character with broken spin degeneracy. A few exceptions 



For Co and Fe in the ZZD structure more bands of 
one type of spin cross the Fermi level as compared to 
those of the other type of spin resulting in a high spin- 
polarization at the Fermi level. This situation implies 
that in the ballistic electron transport, the conductance 
of electrons with one type of spin is superior to elec- 
trons with the opposite type of spin; namely <7j ^> <j\. 
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FIG. 7: (Color online) Energy band structures of 3d-TM 
atomic chains in their L and LD structures. The zero of energy 
is set at the Fermi level. Gray and black lines are minority and 
majority bands, respectively. In the antiferromagnetic state 
majority and minority bands coincide. Energy gaps between 
the valence and the conduction bands are shaded. 



Accordingly, the conductance of electrons across the Fe 
and Co chains becomes strongly dependent on their spin- 
directions. This behavior of the infinite periodic Fe or 
Co chain is expected to be unaltered to some extend 
for long, but finite chains and can be utilized as a spin- 
dependent electronic device. In closing this section, wc 
emphasized that the infinite, periodic chains of 3d-TM 
atoms can be in the zigzag structure corresponding to a 
local minimum. However, most of the zigzag structures 
are dimerized. Dimerization causes remarkable changes 
in electronic and magnetic properties. 



IV. SHORT CHAIN STRUCTURES 

Periodic infinite chains in Section III are only ideal 
structures; long finite-size segments perhaps can attain 
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FIG. 8: (Color online) Energy band structures of 3d-TM 
atomic chains in their zigzag (ZZ) and dimerized zigzag (ZZD) 
structures. The zero of energy is set at the Fermi level. Gray 
and black lines are minority and majority spin bands, respec- 
tively. The gray and dark lines coincide in the antiferromag- 
netic state. Only dark lines describe the bands of nonmag- 
netic state. The energy gap between the valence and the 
conduction bands is shaded. 



their physical properties revealed above. On the other 
hand, the end effects can be crucial for short segments 
consisting of few atoms which may be important for var- 
ious spintronic applications. In this section, we examine 
short segments of 3d-TM chains consisting of n atoms, 
where n=2-7. 



A. Collinear Approximation 

We first study the atomic structure and magnetic prop- 
erties of the finite chains within collinear approximation 
using ultra-soft pseudopotentials4^ The linear structure 
is unstable for the finite size segments. Various planar 
zigzag structures, which are only a local minima, are de- 
scribed in Fig. [TJ We optimized the geometry of these 
zigzag structures with different initial conditions of mag- 
netic moment on the atoms within collinear approxima- 
tion. If the final optimized structures for q different ini- 
tial conditions result in different average cohesive energy 
(or different total energy), they may actually trapped 
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in different local minima. Here we considered following 
different initial conditions: (1) At the beginning, oppo- 
site magnetic moments, ±/x a , have been assigned to adja- 
cent atoms, and the total magnetic moment, ji = X^Mai 
has been forced to vanish at the end of optimization for 
n = 2 — 7. Initial magnetic moment, fi a , on atoms are 
determined from the Hund's rule. (2) For n — 2 — 7, 
initial magnetic moment of all atoms have been taken in 
the same direction, but the final magnetic moment of the 
structure has been determined after optimization with- 
out any constraint. (3) For n = 2 — 7, the system is 
relaxed using spin-unpolarized GGA. (4) For n = 2 — 7, 
initial magnetic moment of chain atoms have been as- 
signed as is done in (1), but \i = ^2 Ha is not forced to 
vanish in the course of relaxation. (5) For n = 2 — 7, 
spin-polarized GGA calculations have been carried out 
without assigning any initial magnetic moment. (6) We 
have assigned the magnetic moment fj, J.| for n = 4, 
and tl I IT for n = 5. Here different spacings be- 
tween two spin-arrows indicate different bond lengths. 
This way different exchange coupling for different bond 



lengths and hence dimerization is accounted. (7) The ini- 
tial magnetic moment on atoms f J. |j f | for n = 6 and 
TIT TIT for n = 7 have been assigned. (8) Similar to 
(7), initial magnetic moment TIT TIT and Tl ITI IT 
have been assigned for n=6 and n=7, respectively. Last 
three initial conditions are taken into consideration due 
to the fact that different bond lengths of 3d-TM atoms 
affect the type of magnetic coupling between consecu- 
tive atoms^ The initial atomic structures have been op- 
timized for these initial conditions except Cu and Zn. 
Only first three conditions are consistent with Cu and 
Zn. As the initial geometry, a segment of n atoms has 
been extracted from the optimized infinite zigzag chain 
and placed in a supercell, where the interatomic distance 
between adjacent chains was greater than 10 A for all 
atoms. Our results are summarized in Table IIII1 where 
the magnetic orders having the same lowest total energy 
occurred p times from q different initial conditions, are 
presented. In this respect the magnetic ordering in Ta- 
ble IIIII may be a potential candidate for the magnetic 
ground state. 



The average cohesive energy of finite-size chains in- 
crease with increasing n. In Fig. [9l we plot the aver- 
age cohesive energy of these small segments consisting 
of n atoms. For the sake of comparison, we presented 
the plots for the linear and zigzag structures. The aver- 
age values of cohesive energy in Fig. [Hb) are taken from 
Table IIIII We note three important conclusions drawn 
from these plots, (i) The cohesive energies of the zigzag 
structures are consistently larger than those of the linear 
structures, and the cohesive energies also increase with 
increasing n. (ii) For each types of structures, as well as 
for each n, the variation of E c with respect to the number 
of 3d electrons in the outer shell, Nd, exhibits a double 
hump shape, which is typical of the bulk and the infinite 
chain structures as presented in Fig. [5] (iii) For specific 
cases E c (n2) < E c (ni), even if n-i > nx (V and Cr). This 
situation occurs because energy cannot be lowered in the 
absence of dimerization. 



Most of the finite zigzag chain structure of 3d-TM 
atoms have a FM lowest energy state. The magnetic 
ordering specified by AFM* for specific chains indicates 
that the magnetic moment on individual atoms, fia, may 
be in opposite directions or may have unequal magni- 
tudes, but the total magnetic moment, fi = ^2/J> a , adds 
up to zero. The finite chains of Zn atom are always non- 
magnetic for all n. Finite zigzag chains of Cu are non- 
magnetic for even n, except n = 6. Interestingly, the 
dimerized linear chain of Cr (n — 5) with a FM low- 
est energy state is more energetic than that of the zigzag 
chain given in Table HTT1 The geometry of this structure is 
such that two dimers consisting of two atoms are formed 



at both ends of the linear structure and a single atom at 
the middle is located equidistant from both dimers. The 
distance from the middle atom to any of the dimers is 
approximately twice of the distance between the atoms 
in the dimer. Even though the nearest neighbor distance 
of the middle atom to dimers is long, there is a bonding 
between them. The cohesive energy is ~ 0.2 eV higher 
than that of the zigzag case, and the total magnetic mo- 
ment of the structure (6 hb) is provided by the atom 
at the middle. This is due to the fact that two dimers 
at both ends are coupled in the AFM order. This is an 
expected result because the cohesive energy (per atom) 
of Cr2 is higher than that of Crs in the zigzag structure. 
The LUMO/HOMO gap for majority and minority spin 
states usually decrease with increasing n. However, de- 
pending of the type of TM atom, the maximum value 
of the gap occurs for different number of atoms n. The 
zigzag chain of Zn atoms usually have the largest gap for 
a given n. 

Even though the total magnetic moment, /i — ^ /i a , of 
the AFM* state vanishes for the finite molecule, LUMO- 
HOMO gaps for majority and minority states are not 
generally the same as in the AFM state. This can be 
explained by examining the magnetic moment on every 
individual atom and the geometry of the molecule. For 
Cr4, the magnetic moment on each atom are lined up as 
described in the sixth initial condition. In this ordering, 
two dimers each consisting of two atoms are in the AFM 
ordering within themselves, but in the FM ordering with 
each other. The distribution of final magnetic moment on 
atoms for Mn6 also obey one of the initial conditions (case 
7) . Three dimers each consisting of two atoms coupled in 
the AFM order within themselves, but in the FM order 
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TABLE III: The average cohesive energy E c (in eV/atom); the net magnetic moment fi, (in Bohr magneton /j,b); magnetic 
ordering (MO); LUMO-HOMO gap of majority /minority states, Eq and Eq, respectively (in eV) for lowest energy zigzag 
structures, p/q indicates that the same optimized structure occurred p times starting from q different initial conditions. (See 
text) Results have been obtained by carrying out structure optimization within collinear approximation using the ultra-soft 
pseudop ot ent ials . 
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0.0 
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with each other. Similar results are also obtained for 
other AFM* states. 

The zigzag planar structure for n > 3 in Table IIIII 
corresponds to a local minimum. To see whether the 
planar zigzag structures are stable, or else it transforms 
to other geometry by itself is a critical issue. To assure 
that the finite chain structures of n = 4 and n = 7 in 
Table Hill are stable in a local minimum, we first displaced 
the atoms out of planes, then we optimized the structure. 

I 



Upon relaxation all displaced atoms returned to their 
equilibrium position on the plane. 



B. Noncollinear approximation and the spin-orbit 
interaction 



In cases where both AFM and FM couplings occur and 
compete with each other, collinear magnetism fails for 
modelling the ground state magnetic ordering. A midway 



between AFM and FM exchange interactions results in 
allowing the spin quantization axis to vary in every site of 
the structure. Geometric structure also influences non- 
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TABLE IV: The average cohesive energy, i? c (in eV/atom); the components (fi x , fi y , fi z ) and the magnitude of the net magnetic 
moment n (in /j.b); LUMO-HOMO gap Eg (in eV) / energy gap under 1 /j,b applied magnetic field along x-direction(/2- 
direction) E G /E G ; magnetic ordering MO; Spin-orbit coupling energy AE x a (/ AE z a ) (in meV) under applied magnetic field 
along £-direction(/z-direction). p/q indicates that the same optimized structure occurred p times starting from q different initial 
contidions. Results have been obtained by carrying out structure optimization calculations within noncollinear approximation 
using PAW potentials. Mn7 is not stable in the planar ZZ structure. 
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FIG. 9: (Color online) Variation of the average cohesive en- 
ergy of small segments of chains consisting of n atoms, (a) 
The linear chains; (b) the zigzag chains . In the plot, the 
lowest energy configurations for each case obtained by opti- 
mization from different initial conditions. 



collinear magnetism. Frustrated antiferromagnets hav- 
ing triangular lattice structure, disordered systems, bro- 
ken symmetry on the surface will result in noncollinear 
magnetism. cv-Mn, spin glasses, domain walls, Fe clus- 
ters are examples of this type. Coupling the magnetic 
moment to the crystal structure (spin-orbit coupling) 
poses the magnetic anisotropy which again results in the 
noncollinear magnetism in the structure. Finite struc- 
tures that is studied in this paper, all have low symme- 
try and AFM-FM coupling competition which increase 
the probability of observing noncollinear magnetism. 
There are many different approaches for implement- 
ing noncollinear magnetism, such as ASW (Augmented 
Spherical Wave), CPA (Coherent-Potential Approxima- 
tion), LSD A (Local-Spin-Dcnsity Approximation). In 
our study, we use fully unconstrained approach to non- 
collinear magnetism. 4 9 The Hamiltonian of the system 
after making simplifications and approximations will be: 



eff T 

J2\Pi>(D?f+ 1 D?f-iD?f)<p j \ 

Here is the effective one-electron potential which 
depends on the electron density and the magnetization at 
each site; — \5 a p stands for the kinetic energy of the sys- 
tem. D"f, 1 D"- i and x Dff terms in the summation sign 
over augmented channels represent the correction terms 
for long range, effective potential and wave functions. For 
further details see Ref . gHHHlsllll 

The finite chains discussed in the previous section 
within collinear approximation will now be treated using 
noncollinear approximation. To this end, the structure of 
chains have been optimized starting from the same initial 
geometry (starting from a segment of n atoms extracted 
from the optimized infinite ZZ (or ZZD) chain placed in 
a supercell) and following five different initial configu- 
rations of spins on individual atoms, (i) The direction 
of the initial magnetic moment on the atoms are con- 
secutively altered as xyzy. (ii) No preset directions are 
assigned to the individual atoms, they are determined in 
the course of structure optimization using noncollinear 
approximation, (iii) For each triangle, the initial mag- 
netic moment on the atoms have a non-zero component 
only in the xy-plane, but (^ fJ- a )xy = 0. (Here 'A' stands 

A 

for the summation over the atoms forming a triangle) 
(iv) Similar to (iii), but (%2fi a ) z / 0. (v) In a zigzag 

A 

chain, the magnetic moment of atoms on the down row 
are directed along z-axis, while those on the up row are 
directed in the opposite direction. Using these five dif- 
ferent initial conditions on the magnetic moment of in- 
dividual atoms, the initial atomic structure is optimized 
using both ultra-soft^ pseudopotentials and PAW— po- 
tentials. We first discuss the results obtained by using 
ultra-soft pseudopotentials. Almost all of the total mag- 
netic moment and the cohesive energy of the optimized 
structures have been in good agreement with those given 
in Table HTT1 (obtained within collinear approximation). 
However, there are some slight changes for specific finite 
structures. For example, SC7 is found to have magnetic 
moment of 7 /is in collinear approximation. Even though 
one of the initial condition in noncollinear calculations re- 
sulted in the same magnetic moment and energy, there 
is even a more energetic state (0.01 eV lower) with total 
magnetic moment of 9 /is (5.2, 5.2, 5.2). The same situ- 
ation also occurred with PAW potential. T15 has a special 
magnetic moment distribution which is the same for both 
ultra-soft and PAW cases and will be explained below. In 
collinear approximation, V5 is noted to have zero mag- 
netic moment, nevertheless there is a state 0.03 eV lower 
in energy which is FM with fi = 1 fi B (0.7, 0.7, 0.2). 
Even though C07 has the same total magnetic moment 
in both collinear and noncollinear case, there is a signif- 
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icant energy difference between two cases. Noncollinear 
structure of C07 has ~ 0.4 eV lower energy with mag- 
netic moment distribution as (8.5, 8.5, 9) [Ib- SC4, SC7, 
Ti 3 , Ti 7 , V 3 , V 7 , Cr 7 , all structures of Mn, Fe 6 and Ni 5 
have truly three dimensional (3D) magnetic moment dis- 
tribution. 

Finally, noncollinear (NC) calculations have been per- 
formed using PAW potentials starting with five differ- 
ent initial assignment of magnetic moments as described 
above. Most of our calculations have yielded the same 
magnetic moment distribution with previous calcula- 
tions, but there are still few cases, which are resulted dif- 
ferently. Mn 7 is an exception; all structure optimization 
starting from different initial conditions resulted in a non- 
planar geometry. Note that in collinear and noncollinear 
calculations using ultra-soft pseudopotential Mn 7 was 
stable in a local minimum corresponding to the planar 
zigzag geometry, but it formed a cluster when spin-orbit 
coupling and NC effects are taken into account. Unlike 
other n — 5 zigzag structures, T'15 has a unique ordering 
of the atomic magnetic moments. Two Ti atoms of the 
upper row have magnetic moments which are in opposite 
direction. Similarly, two Ti atoms at the ends of the lower 
row also have atomic magnetic moments in opposite di- 
rection, but the magnitude of moments are smaller than 
those of the upper row. The atom at the middle of the 
lower row has no magnetic moment. In n = 6 case, only 
C06 has a non- vanishing magnetic moment. Other atoms 
form dimers which are coupled in the AFM order. If we 
assume that the shape of n = 6 molecule is parallelogram, 
there is an AFM coupling between the atoms on both 
diagonals. In addition to these, remaining two atoms in 
the middle also coupled in the AFM order. Atoms in Sc6, 
Tig and Vq have magnetic moments only in the xy-pl&ne, 
whereas Cr@, Mn@ and Ni@ have 3D magnetic moment 
distribution. Cr„ chains exhibit an even-odd disparity; 
Cr„ has an AFM ordering for even n, but it has a FM 
ordering for odd n. There are also cases where collinear 
and noncollinear calculations with ultra-soft pseudopo- 
tential resulted in an excited state for magnetic moment 
distribution. Although PAW potential calculations found 
the same magnetic ordering with collinear and ultra-soft 
NC cases, there are even more energetic states for Sc6, 
V4, Crs and Mns given in Table [TV] Geometric dimer- 
ization also plays an important role in determining the 
average cohesive energy of the system. Interestingly the 
cohesive energies of V6 and V 7 ; Cr@ and Cr 7 ; N14 and 
Ni5 are not changing with n. It should be denoted that 
Hobbs et. al^- carried out noncollinear calculations with 
the PAW potential on Cr 2 _5 and Fe2_5 finite chain struc- 
tures. Here, our results are in agreement with those of 
Hobbs et. al& 

We calculated the spin-orbit (SO) coupling energies, 
A£^ Q and Ai?f , under a unit magnetic field along 
x- and z-direction, respectively. Here, the optimized 
structure of every initial condition together with the 
calculated magnetic moment on the individual atoms are 
used for the calculation of SO coupling. The optimized 



structures of (TM) n and atomic magnetic moments, /x Q , 
have been determined within noncollinear approximation 
using PAW potentials. The results are given in Table llVl 
in units of meV. As it can be easily seen that SO 
coupling does not play an important role on the energy 
of the planar finite structure. However, SO coupling be- 
comes crucial when the total magnetic moments, which 
happen to be oriented in different directions owing to 
the different initial conditions, result in the same energy. 
In this case SO coupling energy difference changes with 
respect to the direction of the applied magnetic field. In 
Table HVl only the most energetic configurations includ- 
ing the SO coupling effects are given. AE^ and AD| 
appear to be independent of n except Mn5, C02 and 
C03. It is also observed that when SO coupling is taken 
into account, LUMO-HOMO gap energies decreases as 
in the bulk structures. Only for N14, Cr5, Mns and Ni 7 , 
LUMO-HOMO gap increased due to the fact that the 
final geometry of SO coupling calculations has further 
relaxed slightly from that of NC calculations. We close 
this section by noting that PAW potential is found to 
be more suitable for the following reasons. First, the 
individual atomic character in the chain structure, as 
well as local magnetic moment are better represented by 
the PAW potential. Secondly, using PAW potential one 
can provide an accurate prediction of spin-orbit coupling. 

V. CONCLUSION 

In this paper, we presented an extensive study of 
the structural, electronic and magnetic properties of 
monatomic chains of light transition metal atoms (Sc, 
Ti, V, Cr, Mn, Fe, Co, Ni as well as Cu and Zn) using 
first-principle plane wave methods. We considered infi- 
nite and periodic chains (with linear, dimerized linear, 
zigzag and dimerized zigzag geometry) and small chains 
including 2-7 atoms. Because of end effects, we found 
a dramatic differences between infinite chains and finite 
ones. Therefore we believe that the basic understanding 
of monatomic TM chains have to comprise both infinite 
and finite structures as done in the present paper. 

The infinite, dimerized linear structures have a shal- 
low minimum only for a few TM atoms; planar zigzag 
and dimerized zigzag structures, however, correspond to 
a lower binding energy providing stability in this geome- 
try. As for short chains consisting of 4-7 TM atoms, the 
planar zigzag structure is only a local minimum. The fi- 
nite chains tend to form clusters if they overcome energy 
barriers. We found close correlation between magnetic 
state and geometry of chain structure. In this study, we 
presented the variation of binding energy as a function 
of lattice constant for different structures and magnetic 
states. We also revealed the dependence of electronic and 
magnetic properties on the atomic structures of chains. 
We found that the geometric structure influence strongly 
the electronic and magnetic properties of the chains. For 
example, infinite linear Vanadium chain becomes half- 
metallic upon dimerization. Similarly, infinite dimerized 



12 



and metallic Sc chain becomes half-metallic with 100% 
spin-polarization at the Fermi level upon transformation 
to zigzag structure. Furthermore, while the infinite linear 
Mn chain has an AFM ground state, with fi = fi a = 0, 
but | Yl Mai = I S Val — 4-40 fiB, it becomes a FM metal 
with /i = A*a = 4.36 /ib as a result of the structural 
transformation from linear to dimerized zigzag structure. 

Magnetic ordering of finite-chains becomes more com- 
plex and requires a treatment using noncollinear approx- 
imation. The structure optimizations carried out using 
ultra-soft pseudopotentials generally result in the same 
cohesive energy and magnetic moment in both collincar 
and noncollinear approximations. However, for specific 
finite chains the total magnetic moments calculated by 
using PAW potentials with the same initial magnetic mo- 



ment distribution differ dramatically from ultra-soft re- 
sults. Of course, our results which covers much more than 
3000 different structure optimizations may not include 
the lowest energy state, but indicates the importance of 
noncollinear treatment. 
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